Measuring billiard eigenfunctions with arbitrary trajectories 
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We propose a method of measuring approximate quantum eigenfunctions in polygonalized billiard 
geometries, based on a quasiclassical evolution operator having a (smoothened) Perron-Frobenius 
kernel modulated by a phase arising from quantum considerations. We show that quasiclassical 
evolution differs from semiclassical (or quantum) evolution but the operators have the same set of 
eigenfunctions. We demonstrate this by determining the quasiclassical eigenfunctions of the polyg- 
onalized stadium billiard using arbitrary trajectories and comparing this with the exact quantum 
stadium eigenfunctions. 

PACS number(s): 05.45.Mt, 03.65. Sq, 31.15.Gy,05.45.Ac 



The quantum billiard problem consists of determin- 
ing the eigenvalues and eigenfunctions of the Helmholtz 
equation V^'0((7) + k'^ijj{q) — with 'ip{q) = on the bil- 
liard boundary dB (Dirichlet boundary condition). This 
simple wave equation arises in various contexts and has 
been used extensively to test ideas of quantum chaos. It 
can describe acoustic waves, modes in microwave cavities 
and has relevance in studies on quantum dots where the 
motion of electrons can be regarded as "free" inside an 
enclosure. The problem is analytically tractable only for 
the small subset of "integrable" boundaries for which the 
classical dynamics is regular. For other enclosures, the 
eigenstates must be computed numerically and a num- 
ber of efficient "boundary" methods exist that allow us 
to study the eigenvalues and eigenfunctions. 

Of particular interest is the determination of approx- 
imate quantum eigenstates using classical quantities. 
While the old quantum theory of Bohr and co-workers 
works only for regular or integrable systems, modern 
semiclassical theories have respondend to the challenge 
posed by chaotic classical dynamics and the successful 
quantization of the Helium atom [Q points to its success. 
The aim however is not necessarily linked to the devel- 
opment of a cheap substitute for the computer intensive 
numerical methods that determine the exact quantum 
states. While this is a desirable consequence, semiclassi- 
cal studies endevour to provide an understanding of the 
quantum phenomenon in terms of classical objects that 
we are so familiar with. Modern semiclassical methods 
have indeed furthered our understanding of the quantum- 
classical correspondence. Thus, we are now aware of the 
duality of quantum eigenenergies and classical periodic 
orbits - a relationship that now forms the cornerstone 
of most semiclassical theories The study of "scars" 
has also revealed the structure of quantum eigenfunc- 
tions and classical trajectories have even been used to 
construct semiclassical eigenfunctions of chaotic systems 

. Quantum states can thus be contemplated in classical 
terms as a first approximation with corrections providing 
its true quantum nature. 

In this context, a question that may be asked is the 
following : what is the degree of classical information 



that is required in order to extract a first approximation 
of a quantum state ? This is especially pertinent when 
the system in question is chaotic or mixed with islands 
of regularity interspersed in the chaotic sea. Since the 
quantum state (or the quasiprobability distributions con- 
structed out of them) can essentially resolve phase space 
structures of the size of a Plank cell, is information finer 
than that redundant ? For billiards, the de Broglie wave- 
length A, provides a relevant length scale that can be used 
effectively to probe the boundary of the enclosure. If a 
smooth billiard boundary is polygonalized such that the 
short time classical dynamics is well approximated, the 
two billiards are semiclassically equivalent, provided, A 
is larger than the average length of edges of the poly- 
gon 01 . Thus, instead of the full chaotic dynamics of the 
stadium billiard, one may as well consider the dynamics 
of its polygonal counterpart for a given de Broglie wave- 
length. The polygonalization concept is implicit even 
in the numerical computation of the "exact" eigenvalues 
using boundary methods in which the perimeter is dis- 
cretized with the number of points N ^ C/X where C is 
the perimeter of the enclosure. This idea has however not 
been used to determine semiclassical eigenvalues or eigen- 
functions. Indeed, by most accounts a polygonalized ap- 
proach to semiclassics is bound to be even more difficult 
since periodic orbit quantization of polygons has proved 
largely unsuccessful Also, it is generally believed that 
diffractive contributions must be included even for ob- 
taining a first approximation of a polygonal quantum 
state. On the other hand, it has recently been shown 
that closed alomost-periodic (CAP) orbits also con- 
tribute in generic polygonal enclosures with weights that 
are comparable to periodic orbits. Besides, they are more 
numerous and hence indispensable for semiclassical quan- 
tization. Thus, the failure of periodic orbit quantization 
in polygons is not so much due to diffraction as due to 
the neglect of CAP orbits. The modified periodic orbit 
theory for polygons should in principle be able to provide 
approximate quantum eigenvalues although CAP contri- 
butions can be difficult to evaluate. 

A recently developed time domain technique of 
determining quantum eigenvalues in marginally stable 
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billiard geometries makes use of arbitrary classical tra- 
jectories. The algorithm involves shooting trajectories 
in various directions from a point interior to the billiard 
(call it q') and at each time step, recording the fraction 
of trajectories that are in an e neighbourhood of a point 
q, weighted by a phase arising from quantum considera- 
tions. The peak positions in the power spectrum of this 
weighted fraction, F[t)^ are related to the quantum eigen- 
values and as we shall show in this letter, the heights of 
the peaks are a measure of the quantum eigenfunctions at 
the point q. Thus, the arbitrary trajectory quantization 
method (ATQM) is aptly suited for polygonalized bil- 
liards and we shall demonstrate for the stadium billiard 
that the quasiclassical eigenfunctions of the polygonal- 
ized stadium do approximate the "exact" eigenfunctions 
of the full stadium. 

The ATQM bypasses the need to evaluate periodic 
or CAP orbits but is nevertheless based on the modi- 
fied periodic orbit theory for polygons. The parameter 
e (©(l/fc)) embodies the contribition of CAP orbits and 
must be nonzero for any calculation at finite k. In the 
following, we shall first use a plane wave ansatz to show 
that the peak heights in the power spectrum are a mea- 
sure of the approximate quantum eigenfunction and then 
demonstrate this for a polygonalized stadium. 

We define the quasiclassical |9j propagator C*^^{(p) re- 
stricted to the ip invariant surface as 
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dq'%-g'*((p))e— 0(g') (1) 



where g'* (ip) is the position at time f of a trajectory which 
starts at q' {t = 0) on an invariant surface labelled by (p 
and the energy E. The phase ^{t) = v{q'^{if)) depends 
on the caustic structure of the trajectory and is identical 
to the phase in the semiclassical propagator ||^. When 
v{t) is identially zero (as in case of Neumann bound- 
ary conditions), C*q^{if) reduces to the Perron- Frobenius 
operator on the invariant surface. Note that rational 
polygonal billiards have a second constant of motion due 
to which the invariant surface is 2-dimensional 
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successive reflections, the trajectory changes its momen- 
tum following the laws of reflection but continues to live 
on the same invariant surface labelled by E and p. Thus, 
Lp also labels distinct trajectories that start from a given 
point q' . 

The full quasiclassical evolution operator is defined as 



dq'l [ dp, %-g'*(^))e--(*W2^0(g') 

(2) 



= J dq Kqciq,q',t) (j){q) 



Note that quasiclassical evolution differs from semiclas- 
sical evolution. To illustrate this, consider a particle in 
a one-dimensional box (Dirichlet boundary conditions) 
and consider the evolution of the quantum eigenfunction 



ipn{q) = e*'^"' — e~*'^"*, kn = nir/L . Its time evolution 
in quantum mechanics is simply e~*^"*'0(9) where En = 
h'^kf^/2m. Its quasiclassical evolution, £g^(-f ) o -tpniq), is 
given by 



(3) 



where n{q^*{+v)) is the number of reflections suffered by 
a trajectory in time —t with initial position q and initial 
velocity +v. Similarly, £*^(— ) o ip^^q) is 



-ik„{q '(—"))! p-i-^n{q '(-f)) 



(4) 



with the — sign in £*^(— ) denoting negative velocity. 
Note that the flow is such that the velocity changes sign 
at every reflection from the walls at q = and q = L while 
n{t) increments by one at each of these instants. For the 
flow q~'^{+v), the reflections occur at = {q + nL)/v 
so that for < t < tf, q~\+v) = v{t ~ t^) = vt - q. 
Similarly, for the flow (j~*(— w), the reflections occur at 
t:^ = (L — q + nL)/v and for <q < t < , q~^{—v) — 
L — v{t — Iq) = 2L — vt — q. It follows hence that 



£* (±)oV^„(q) = e^ 



for all t. Thus 



-ikn{q =p vt) 



(5) 



(/:*,(-f) -f £*,(-)) o^„(g) = 2cosiknvt)Mq) (6) 

In other words, ipniq) is also an eigenfunction of the full 
quasiclassical evolution operator = £*^(+) + £^^(— ). 
Note that in the Neumann case, ipnil) is not an eigen- 
function; rather e'*^""^ -I- e~*'^"^ is an eigenfunction with 
n{t) = 0. 

For a general 2-dimensional billiard, there is strong 
numerical evidence to suggest that a (real) plane wave 
superposition can be used to construct eigenfunctions 
pTj-p^ . We shall adopt here the view that a flnite plane 
wave superposition does yield at least good approximate 
eigenfunctions. For polygonalized billiards, the semiclas- 
sical wavefunction can be expressed as 
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E 



ik cos( fij)x-\-ik sin[fij)y 



(7) 



where Aj are constants ]14| and the number of terms M 
in the expansion is determined by closure of the wave 
vector k — (fc cos /ij, fc sin/ij) under reflection from the 
edges. 

For this finite superposition of plane waves, the bound- 
ary condition ip{q) =0 on dB can be satisfied if the waves 
vanish in pairs with an incident wave giving rise to a re- 
flected wave. Thus on the Ith segment y = aix + bi, we 
must have 



A ^i(k cos +a; k sin fij)x-\-ibi k sin 
_^ ^ ^i(k cos +aik s'm fi,j,)x+ibik sin fiji _ q 



(8) 
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Assuming that Hj' is related to jij through the laws of 
reflection, it is easy to show that 

cos /ij + ai sin fij — cos /ij' + a; sin /i^/ (9) 

Thus, eq. (||) reduces to 

where /ij' ~ n — fij + 29i and 0; is the angle between the 
positive X-axis and the outward normal to the ^th line 
segment. 

Note that for each of the K segments on the boundary, 
the jth wave has in general a different reflected wave as a 
counterpart so that eq. (|l^) gives K different expressions 
for Aj . In general (barring exceptions such as the rectan- 
gle billiard), these "boundary conditions" can be satisfied 
only approximately as we shall argue below. Recall that 
for the numerical determination of exact eigenvalues us- 
ing a plane wave basis, the boundary is discretized {N 
points) and an appropriate measure (such as a determi- 
nant) is used to determine the eigenstates which satisfy 
the boundary condition at these points. Convergence 
can be achieved by increasing N so that as TV — > oo, the 
boundary condition is satisfied exactly. In contrast, the 
number of terms in eq. |^ is fixed. Thus, if the exact eigen- 
function contains additional plane waves, the boundary 
condition will be satisfied approximately and the plane 
wave expansion of eq. |^ can only give an approximate 
quantum eigenfunction. 

We shall now establish that the (finite) plane wave su- 
perposition (eq. 1^) is also an approximate eigenfunction 
of the quasiclassical evolution operator provided the set 
of "quantization" conditions given by eq. ( p^ are satis- 
fied. 

Consider therefore the plane wave superposition of 
eq. (^. Its quasiclassical evolution is given by 
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'(¥>)+ifcyi/ *(ip)g-m(t) 



(11) 



where = A:cos(/ij), ky = ksm{fij) while x~*(ip) and 
y~*(v) denote the fiow at time — t with initial position 
{x,y) and velocity {v cos (p,v sin ip). For short times, this 
is given by 

M 

C} ((^) = yi^.g»fex(2:--ucosipt)+*fey(i'-i'sini^t) ^y^^ 

J = l 

As before, we shall first determine the evolution of a sin- 
gle wave after reflection from one of the segments, y — 
ai.x + hi.. For the flow, (a;^*((^), ((/?)), reflection from 
the line segment takes place at — {x — a;o)/(f cos(/?) — 
{y ~ yo)/('i'sin(/3) where (a;o,j/o) is the point of impact. 
The flow at a time t after the reflection is given by 



X *{(p) = x{t) ^ xq + V cos{lp — 26i){t ~ to) 
y-\ip)^y{t)^yo + vsin{Lp-26i){t^to). 



It is easy to verify that after one reflection from the seg- 
ment y = aix + bi, the wave Aje^'^^"'^'^^^'^''^^'^'^^^ evolves 
quasiclassically to 

^ ^ cos fj,^/ {x~v cos :fit)+k sin /j^/ {y — v sin ^pt) (14) 

where /ij' = tt — /ij +20/ and 

^^.gifcf,, sinM, ^^,e*fe6i ^"^1^,' ^ (15) 

Thus, after one reflection, the finite plane wave superpo- 
sition assumes the form for small t provided the reflected 
waves are included in the superposition and the "quan- 
tization conditions" are (approximately) satisfied for a 
given value of fc„. It follows that '4>n{<l) is an (approxi- 
mate) eigenfunction of Cqc = J dip Cqc{ip) : 



dip Cqc{p)'4'n{q) 



(16) 



= 2TTJo{knVt) Ml) (18) 

where Sj{kn) — kn cos{fij)x + kn sin(/ij)y. Thus, a finite 
plane wave superposition can be an approximate semi- 
classical and a quasiclassical eigenfunction under identi- 
cal conditions. 

We shall now demonstrate our result for a stadium bil- 
liard consisting of two parallel straight segments of length 
2 joined on either end by a semicircle of unit radius. 
For the evaluation of the quasiclassical eigenfunctions, we 
shall consider a polygonalized enclosure where each semi- 
circle is replaced by 12 straight edges of equal length. In 
order to determine the quasiclassical eigenfunctions, we 
shall first evaluate a smoothened quasiclassical kernel 



Kqc{q,q',t) 



dip Je('7-g'*(<^))e~""(«"('^)) 



(13) 



= J2Mq)i'n*iq')Anit) (19) 

n 

as a function of time. Thus is achieved by shooting trajec- 
tories from a point q' at various angles and evaluating the 
fraction of trajectories in a cell of size e jl^ at q, weighted 
by the phase e-*'^"(«"('^)). Since A„ = 27rJo(fc„wi), for 
ti = 1, a fourier transform of Kqc{q, q' ,t) has peaks at 
k = kn and the heights are proportional to ipniq)- 

Note that the smoothening of the delta function kernel 
is essential in order to accomodate closed almost-periodic 
orbits and shows up naturally in an alternate proof (of 
the identity of the semiclassical and quasiclassical eigen- 
values) involving the trace of the quasiclassical and semi- 
classical propagators In the present formalism in- 
volving plane waves, smearing of the kernel leads to a 
modified "quantization condition" for the quasiclassical 
eigenfunctions. We shall however ignore these complica- 
tions and merely reiterate that the parameter e ~ 1/k. 

Fig |l|a shows a "bouncing ball" quasiclassical eigen- 
hmction intensity |V'n(9)P at fc = 10.97 in the quarter 
stadium. 
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FIG. 1. (a) A quasiclassical bouncing ball eigenfunction of 
the polygonalized stadium and (b) its quantum counterpart 
in the smooth stadium 
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FIG. 2. (a) A quasiclassical eigenfunction peaked along the 
X axis and (b) its quantum counterpart 



stadium at = 11.05 is shown in figj^b. An example of a 
quasiclassical eigenfunction at A: = 4.02 peaked along the 
X-axis is shown in fig. ^ along with its quantum conter- 
part at fc = 4.38. The quasiclassical eigenfunction clearly 
provides a first approximation of the quantum eigenfunc- 
tion in both cases. Eigenfunctions of other billiards in- 
cluding triangles have also been obtained. Details of this 
work will be published elsewhere. 

In conclusion, we have demonstrated that the quasi- 
classical eigenfunctions determined using arbitrary tra- 
jectories in a polygonalized chaotic enclosure, approxi- 
mates the quantum eigenfunctions of the smooth billiard. 
We have also shown that a finite plane wave expansion is 
an approximate eigenfunction of the quantum and quasi- 
classical evolution operators under identical conditions. 
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